home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / cone.f < prev    next >
Text File  |  2002-06-11  |  691b  |  29 lines
  1. C
  2. C PAGE 348: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: CONE.FOR
  5. C
  6. C ICE CREAM CONE EXAMPLE (RANDOM)
  7. C
  8.       DATA  L,M,N/256,0,1000/ 
  9.       DO 2 K = 1,N
  10.         X = 2.0*RANDOM(L) - 1.0 
  11.         Y = 2.0*RANDOM(L) - 1.0 
  12.         Z = 2.0*RANDOM(L) 
  13.         IF( (X*X + Y*Y .LE. Z*Z) .AND.
  14.      A      (X*X + Y*Y .LE. Z*(2.0 - Z)) )  M = M + 1     
  15.         IF( MOD(K,100) .EQ. 0 )  THEN 
  16.           VOL = 8.0*REAL(M)/REAL(N) 
  17.           PRINT *,VOL 
  18.         END IF    
  19.    2  CONTINUE    
  20.    3  FORMAT(5X,I5,F10.5) 
  21.       STOP
  22.       END 
  23.   
  24.       FUNCTION RANDOM(L)
  25.       L = MOD(16807*L,2147483647)
  26.       RANDOM = REAL(L)*4.6566128752458E-10      
  27.       RETURN
  28.       END 
  29.